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Abstract 

The factorial moments (FM) of multiplicity distributions in one- and 
two-dimensions are studied for Pb+Pb collisions at 158 AGeV/c. The 
results are compared with FPJTIOF, RQMD and VENUS predictions. In 
the relation between FM and the number of partitions, M, obvious pos- 
itive intermittence exponents, <f> q , are observed. The exponents, cj> q , are 
compared with the anomalous dimensions, d q . Then they are expressed in 
terms of Levy stable indices, fi. Meanwhile the experimental data within 
the interval, 2 < n < 6, fulfills the requirements of the Levy stable region 
(0 < ii < 2), the results from the other intervals are away from this region. 
All their corresponding indices have negative values. 

PACS: 25.75. G, 61.43, 12.38.M, 64.60. Q 



1 Introduction 

Since the observations of JACEE collaboration Q and the pioneer works of 
Bailas and Peschanski |2|, |3j around mid of 1980's, the intermittent behavior 
in multi-particle production has increasingly gained full credibility among the 
high-energy physicists. The concept "intermittence" itself is primarily used to 
describe the fluctuations observed in the particle density distributions. In many 
experiments, the large dynamic fluctuations at the final state of particle produc- 
tion has been confirmed. To get a satisfactory explanation for the power-law 
behavior of FM with the partition number, M, many theoretical approaches 
are proposed. For example, self-similar branching, QCD parton cascade, multi- 
particle cascade, second-order phase-transition, etc. Also the effects of Bose- 
Einstein correlations on the two-particle correlation functions ^, ||, ^ |7| may 
explain this behavior. The randomization and the hadronic Cerenkov radiation 
are suggested to be partially responsible for such power-scaling behavior. Also, 
if QGP has been produced, the effects of the soft and hard collisions, the cluster- 
ing and resonance decay and the showering in the QGP-surrounding manifest 
themselves in forms of: multiplicity fluctuation, power-scaling behavior, self- 
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similar branching, jet structure, etc. Some of these effects are responsible for a 
kind of canonical fluctuations in the particle production. 

In this paper, the behavior of FM is studied for Pb+Pb collisions at 158 AGeV/c 
with successive partitions in one- (pseudo-rapidity 77 or azimuthal angles ip £ 
{0, 27r}), and two-dimensions (77 and </>). The produced particles are classified 
into four groups according to their pseudo-rapidity values (3 < 77 < 4, 3 < 77 < 5, 
2 < 77 < 5 and 2 < 77 < 6). The observed intermittence exponents, <p q , are com- 
pared with the anomalous dimensions, d q (Renyi dimensions, TZ q ). After that, 
they are plotted in the Levy space. 

This article is organized as the following: The experimental material and the 
simulation by using the three events generators, FPJTIOF, RQMD and VENUS 
are complemented in Section ||. The definitions and formalism are given in 
Section [|. The results of FM and the intermittence exponents, <p q , in one- and 
in two-dimensions are introduced in Sections ^ and ^, respectively. The relations 
between the intermittence and the anomalous exponents and their expressions 
in terms of Levy indices, fi, are discussed in Section ^. Section ^ contains some 
remarks and the final conclusions. 

2 Experimental Material 

The data used in this paper is retrieved from some of the Pb-chambers irradi- 
ated at CERN-SPS during 1996 for the EMU01 collaboration. This material was 
consummately measured and analyzed at the Marburg University by using our 
measuring system MIRACLE Lab ||, || [n^. The collisions have been recorded 
with 82 Pb-beam accelerated to incident momentum of 158 AGeV/c and then 
fired towards stationary lead foil with thickness ~ 250 fim. The lead foil was 
positioned in the front of seven consecutively arranged plastic sheets coated 
on both sides by the nuclear emulsion, FUJI ET-7B. Because of the relatively 
high transverse momentum, most of, if not all, produced particles and frag- 
mentations are emitted within a narrow forward cone. Therefore, their tracks 
can be registered within the forward emulsion sheets. The exposure process is 
mainly controlled by counting the heavily ionizing particles by using scintilla- 
tor and discriminator with relatively high threshold settings. The number of 
beam-particles is registered by using an additional counter installed behind the 
emulsion chamber. Scale and driver electronics (CAEN N145) are used to build 
up pulses to be transferred to the SPS- "kicker" magnet. The kicker magnet con- 
sequently removes the beam, if the particle number transcends 3000 per spot. 
The beam density is ~ 5 • 10 2 nuclei/cm 2 . 

The scanning efficiency in the emulsion chambers is 0.75 ± 0.05 |llj]. The emul- 
sion sensitivity for singly charged particles is as good as 30 grains per 100 
fim. Depending on the incident momentum, the polar angle can be calcu- 
lated. For the actual data, we get polar angle of 1.3 mrad. Depending on 
topology of the microscope's field of view, particles with space angles, 9 < 30° 
(77 = — lntan(0/2) > 1.32) can be acquired. Obviously, most of produced parti- 
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cles are singly charged. It is expected that they are mixed with contamination of 
e~ — e + pairs from Dalitz decay and 7-conversion. The possible over estimation 
of the observed particles density has been determined as ~ 2% || . As a reason of 
the reconstruction algorithm jl(| applied for the measuring system MIRACLE 
Lab, the tracks of these electrons are completely disregarded. In comparison 
with human operators, it is found that MIRACLE Lab has an efficiency up to 
96% Q. Large part of this 4% discrepancy is obviously coming from the par- 
ticles emitted with large space angles (6 > 30°). The missing measurements, 
frequent scattering, unresolved close-pairs and pair production represent addi- 
tional sources for this discrepancy M. Therefore, we could consciously suggest 



to renounce the discussion of the effects of 7-conversion on FM (see Section 4.2 
and Section 5.3 for more details). 



2.1 Simulation of the experimental data 

The investigation of FM is performed for the collisions with restricted particle 
multiplicity^ > 1200. From these events, only the particles emitted within 
predefined pseudo-rapidity intervals (3 < 77 < 4, 3 < 77 < 5, 2 < 77 < 5 and 
2 < ?y < 6) are taken into account. These 77-intervals are chosen at and around 
the mid-rapidity region. Due to the chiral symmetry breaking, the produced 
particles are mainly pions, 7r's. Therefore, the restriction on the ^-values is 
necessary to study the phase-transition, which, in turn, should be restricted to 
the produced particles only. Therefore, the narrow ^-intervals are the suitable 
ones to study the produced particles and consequently the phase-transition. 
Once again, here we used another restriction on the considered phase-space.^ 
On the azimuthal angles, there is no restriction at all, (for each Ar/ group, 
the considered 0-window is allowed to take any real value within the available 
spectrum, {0, 2tt}). 

For the event generators FRITIOF 7.02 ||, RQMD 2.1 @ and VENUS 4.12 
[ fL5| , the input parameters are adjusted for the destination to produce simulated 
events with total multiplicity not smaller than 1200. Each simulated event has to 
go through different cutting processes, in order to simulate, as good as possible, 
the "real" events (Section ||). After these operations, the number of particles 
observed within the space angle, 9 < 30°, has to be > 1200. The criteria used 
for the data cutting, can be taken from |8|, [Tc| ]: 

1. P = p/E > 0.7. In practice, we can mainly measure singly charged par- 
ticles. More than 90% of them are charged 7r's. The number of highly 
charged particles is relatively small jll| . 

2. 77 = — In tan(0/2) > 1.32. The space angle is restricted within 6 < 30° de- 
pending on the geometry of the emulsion chamber, and of the microscope's 
field of view. 

1 This is a restriction on the multiplicity (see Section jj] for the effects of the multiplicity 
restriction on FM [|1J 

2 The influences of the restriction on the considered phase-space to get flat distributions 
of FM are discussed in Section H. 
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3. Charged particles. The emulsion is sensitive for the charged particles only. 



3 Definitions and Formalism 

If the phase-space, Aij, is split into M bins of equal sizes, Srj — Arj/M, the 
exclusive scaled factorial moments [|], [| are given as 

^ n(n - 1) • • • (n - q + 1) 

n m is the multiplicity in m-th bin and n is the average multiplicity in the 
whole A?i-window. With successive partition, At? — > 0, a power-law scaling can 
describe the fluctuations in the particle production, 



q n™ , 

i—1 rjr, 



F q (M) = M*-> £ 'J , (2) 

m=1 II I dTH(&\m,V2,-~,ri q ) 

Pn\- ■ ■) is the pseudo-rapidity density in a sub-interval with q particles investi- 
gated in an event with n total multiplicity. t]q is the minimum pseudo-rapidity in 
m-th bin and rj m is the pseudo-rapidity value in this bin (ry m = T]q + m ■ At]/M). 
Then the factorial moments in m-th bin are to be given as 

q Vm , , 

II / din p)? 1 (r]i, r?2, •••,%) 
FJM) = i=1 Vm ~l . (3) 

In the case that there is no correlation between the produced particles, FM, in 
the whole Ary-interval, can be given as 

F q (M)=M«- 1 Y,[ J dr, m P n {n)\ , (4) 



m=l 



Pn(v) = Pn{ii)/ n is the normalized distribution of pseudo-rapidity density in a 
bin of with size 77. When rj — > A77, we get the maximum value 

F,(M) max = AiT 1 / dr, [P N ( v )] q . (5) 

no 

When F q (M) max — 1, the single-particle distribution in the A77- window will be 
a plane one^| (a distribution with non-statistical fluctuations). 

3 A powerful method to produce flat particle density distributions from non-flat ones can 
be taken from [hql. 
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In this section, we have discussed the exclusive FM and their limitations to 
describe the multiplicity fluctuations. In the next section, we pass to their 
possible applications to describe the particle production in heavy-ion collisions. 



3.1 Factorial moments and intermittence exponents 

For many reasons^, the direct calculation of the exclusive FM is not easy, espe- 
cially for high statistics. Therefore, it has been suggested ||, [| |l8| to implement 
another kind of FM, namely the inclusive FM, 

FJM)™ 1 = A/*- 1 Y Mn m ~l)---(n m -g+l)) ^ (g) 

m— 1 

The value, n, is coming out by dividing the observed multiplicity by the bin 
number, M , and then by the number of events, N. The exclusive and inclusive 
FM are the same, if the particle production can be described by Poisson func- 
tion. In principle, the particle density distribution depends on the dynamics of 
the interacting system (mass, energy, number of events, whether there is any 
restriction in the data sample, 6 etc.). 

pincl n(n-l)---(n-q+l) excl 

Then in the real case, the difference between both types of FM is actually the 
long-range correlations, which appear as a reason of the strong effects of the 
multiplicity on the correlation functions jl2], |l9[ . 
By definition §, § f| FM used in the particle production are 

< F q (M) >= / < nk {Hk ^ [ l k q q + 1} > V (8) 

< nk > is the average multiplicity in fc-th bin of size dn. According to the self- 
similar density fluctuations, the successive partitions of the pseudo-rapidity, Arj, 
leads to the following power-law dependency 

< F q {M) > cx M^ q . (9) 

The exponents, <p q , are called the "intermittence exponents", which practically 
can be retrieved from the slopes of the relations between log F q and log M. They 
are also related to the "anomalous dimensions" |^| 

d q = (10 ) 

The fractal Renyi dimensions |23] are defined as functions of the "anomalous 
dimensions" d a 



4 Mainly, there are two reasons: the limitation on the considered pseudo-rapidity, in order 
to get flat particle density distributions and the limitation on multiplicity JL7| |, in order to 
avoid any possible dynamical correlation between the produced particles. 
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n q = n ■ (i - d q ), (ii) 

The constant, TZ, is the topological Renyi dimension. For multi-fractal processes, 
d q can be given as functions of the orders, q, only (Section |6.3| and Section 
6.4). For mono-fractal processes (Section^), d q are constant. Actually, Renyi 
dimensions can be used to measure the randomization in the produced particles. 
For the total randomization, <p q = 0, d q = 0, and consequently lZ q = TZ. 



4 Results in One-Dimension 
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Figure 1: Tfte factorial moments, F q , for the orders, q = {2,3,4}, are plotted as 
functions of small partition number, M v , in n-dimension. The data used here are 
Pb+Pb collisions at 158 AGeV/c with total multiplicity > 1200. The different Ar/- 
windows are applied above the whole region of pseudo-rapidity (n > 1.32). Only the 
secondary particles emitted within these windows are taken into consider. The data 
are fitted as straight lines. Their slopes are the "intermittence exponents" <j>q- 



In both "real" and "simulated" collisions, pseudo-rapidity intervals are succes- 
sively split into M equal sub-intervals (bins). In each such bins, the multiplicity 
and then the corresponding q-order FM are calculated according to Eq. (|J). 
Fig. |l| shows F q as functions of M v fo r th e orders, q — {2,3,4}, and for the 



different Ayy-intervals given in Section 2.1. The underscore in M v means that 



the partition has been performed in ry-dimension. The splitting in the other 
dimension is completely avoided here. We notice that all pictures show obvious 
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Figure 2: TTie factorial moments, F q , for the orders, q = {2, 3, 4}, are plotted as func- 
tions of small partition number, M^, in ^-dimension. As in Fig. [1] the data are fitted 
as straight lines. Correspondingly their slopes give the component of "intermittence 
exponents" in (^-dimension, (j> q . 



positive intermittence slopes, <p q . Also, it is to realize that the FM- values in- 
crease with the decreasing of A^-intervals. Another remarkable finding noticed 
here is that the slopes increase with the increasing of the orders, q (review Eq. 
(ph). The analysis depicted in Fig. ^ is for the partition in (^-dimension. Nearly, 
the same qualitative features as in previous figure are noticed here. Through 
increasing the partition number, M v or M$, the values of the corresponding FM 
raise fast, as shown in Fig. ||andFig. [|. The reasons for this upwards increasing 
will be discussed later (Section |i~l|). For now, the discussion is restricted only on 
the phenomenological behavior of the increasing of FM with M {vW . In Fig. | 
the smallest 5r] in all pictures is ~ 0.01, which seems to be coincident with the 
detector resolution. Beyond this value, FM reduce suddenly (not shown here), 
i.e. the upwards increasing of FM (shown here) gets inverted. This boundary 
limit can be taken as the low- value of the detector resolution. Until such point, 
we can effectively count the produced particles and then correctly calculate their 
FM. Beyond it, it is impossible to assure whether the measured particles are 
real ones and consequently their calculated FM must be associated with large 
statistical errors. The range of the linear dependence in Fig. || becomes larger 
with wider A?7-windows. Some kind of plateau begin to appear with the in- 
creasing of Arj. Figure ^ depicts the same investigations as Fig. ||, but for the 
partition in the (^-dimension. Nearly, same qualitative features as in Fig. are 
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Figure 3: The factorial moments, F q , for the orders, q = {2,3,4}, are depicted as 
functions of large partition number, M v , in rj-dimension. The partition process is 
repeated as long as the FM-values increase. An exponential growth of FM is noticed 
for relatively large M v . Within the largest Arj-interval, we notice that the increasing 
of FM is saturated. 



noticed here. The only difference is that there is no plateau in any curvature 
region. 

In addition to the large statistical errors, there are other reasons for the negative 
dependence of F ? (M{,ui) on large Mi^y. The properties of FM themselves. 
The FM used here can be applied, only if the following conditions are simulta- 
neously fulfilled: 

• The order q has an integer positive value. 

• The multiplicity in each bin (n m ) is larger than or at least equal to q. 

• q is greater than or at least equal to 2. 

In Fig. U and Fig. |4|, we could expect that the multiplicity, be smaller 

than q. If the limit, q < n m , has been reached, we cannot correctly implement 
Eq. (^). The discussion of FM for large M^uy should not, in the one hand, lead 
to the conclusion that we would like to go through any argumentation of the 
detector resolution. But, as we think, it is a worthwhile finding to realize that 
our measuring system MIRACLE Lab J§|, ||, [l(| is effectively able to determine 
the emulsion resolution. This efficiency has been estimated (for example, by 
testing its reproduction of the measured tracks, when the emulsion plates are 
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Figure 4: The factorial moments, F q , for the orders, q = {2,3,4}, as functions of 
relatively large partition number, M^, in (^-dimension. As in Fig. the partition 
process is repeated as long as the FM-values increase. Here there is no evidence for 
the saturation region. 



randomly rotated). In the other hand, this discussion enabled us to refer to the 
limitations of the used FM. 

The predictions from the three events generators are illustrated in Fig. [5]. For 
FRITIOF simulation one notice that FM have non-linear dependencies for rela- 
tively small M v . Also, the values of FM are not coincident with the experimental 
ones shown in Fig. [|. For large M v , the upwards increasing is to be observed 
until bin size Sr/ — > 0.005. Within these small bins, large numerical uncertainties 
in the computations are obviously expected. Therefore, these regions should be 
charily studied and no serious conclusions according to their features should be 
drawn. The saturation plateau here is too weak. Obviously, they begin later 
and have smaller width than the experimental ones (Fig. ||). From the results of 
RQMD calculations nearly the same features as in previous picture are observed. 
The width of the saturation plateau is larger than that in FRITIOF. Although 
the saturation regions of the VENUS calculations seem to be comparable with 
the corresponding regions in Fig. ^, they are still unable to completely predict 
the experimental results. Obviously, the different discrepancies between the ex- 
perimental data and the three models become larger with increasing M v (see 
Fig. [ll], later). This may be explained due to the effects of Bose-Einstein cor- 
relations, which are not included in any model. In addition to this, during the 
experimental observations p3], these effects have not been considered. Another 
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Figure 5: The FM from FRITIOF, 
RQMD and VENUS simulations 
depicted as functions of M v . The 
left column in each picture depicts 
the partition for relatively small 
partition number, M v . Obviously 
the experimental points cannot be 
fitted as straight lines. The right 
column gives FM for large M v . In 
all pictures an exponential increas- 
ing is noticed. In RQMD and 
VENUS simulations the increase of 
FM is saturated especially in the 
largest Ar/-interval. But the results 
FRITIOF show no evidence for any 
saturation region. 



source for the discrepancies may be the physics of the phase-transition which 
has not been taken into account in any model. 
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4.1 Upwards increasing 

Going back to the gigantic upwards enhancement of FM with large Mjj^, We 
would like to begin here with two points: First, the results published in 25 , in 
which the authors claimed that there is no more increasing in FM in the regions 
beyond the limit dr) < 0.1. In this present work, we clearly notice that FM still 
have an obvious linear increasing, even up to the partition width Srj ss 0.02 (sec 
Fig. and Fig. |). Beyond this narrow bin size, the increasing of FM does not 
only get ahead, but it becomes exponential. Second, a reasonable explanation of 
the tremendous growth of FM for large M^^y may be based on the short-range 
conventional correlations |^6| , on the hadronic resonance and cluster decay j|, || , 
on the Bose-Einstein correlation, etc. In the following is another reason for this 
exponential increase |2?j, ^8|: Using trivial geometric model, each projectile nu- 
cleon interacts along the cylindrical tube in the target nucleus with more than 
one nucleon. Therefore, for each of these sub-nucleus interactions, there is a 



pseudo-rapidity distribution |27 with certain mid-rapidity depending on the 
typical characters of such interacting systems. It depends also on the corre- 
sponding cross section, which in turn depends on whether the knock takes place 
between nucleons in face-to-face or between eclipsed ones [ p9[ . In addition to 
these sub-interactions, the secondary particles (first and/or higher generations) 
produced within the target material cause further interactions with each other 
and with the target and projectile materials. These secondary interactions pro- 
duce additional pseudo-rapidity distributions with their own mid-rapidity cen- 
ters. Therefore, for each of these elementary collisions, a sub pseudo-rapidity 
center can be achieved. The superposition of these elementary pseudo-rapidity 
values taken place in sw&-intervals dr], in which the original intervals Ar/ has 
been partitioned, makes the total effective partition number, M eff , larger than 
M (M cff >M) 0. The reasons for these additional effects can summarize as 
the following: 

1. The successive elementary collisions cannot completely be coincident with 
each others. Therefore the final summation of sub pseudo-rapidity centers 
(X) ^v) produces one distribution with many centers (M eff ^> M). 

2. Even in the ideal case: The superposition produces only one distribution 
with only one center (M off = M), its position has not to be coincident 
with the center of original one, Ar/. In this case, M oS ^ M. 

Evidently, these facts are able to explain the exponential increasing of FM (see 
Fig. H later). The effects of M eff cannot be neglected in such large interacting 
system like Pb+Pb, since the degree of upwards increasing strongly depends on 
the number of participating projectile nucleons and on the size of the target 



tubes | 29 



4.2 Influence of 7 conversion on ID FM 

The measured FM are expected to be partially influenced by the e~-e + — pairs 
produced by Dalitz decay and/or by 7-conversion. In |q], the probabilities of 7- 
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conversion in the emulsion and in the lead foil arc calculated and found that the 
percentage of 7's converting to electron-pairs in the lead foil is ~ 1.71 ± 0.1 and 
in the emulsion is ~ 0.48 ±0.06. Taking the detector resolution and the possible 
space angles between such pairs into account, it is possible to estimate that 
> 78% of these electrons cannot be grabbed by MIRACLE Lab. In addition 
to this effect of unresolved close pairs, most of, if not all, the rest of these 
electrons are obviously included within the so-called 4% discrepancy (Section 
||). Therefore, we could consider that there is a negligible effect of 7's on ID FM. 
Another argument to support this conclusion is coming out from the simulation 
processes. Through adding 7's to FRITIOF for example, it was found f^pf : 
First, the influence from the Dalitz decay is obviously negligible. Second, the 
influence from the 7-conversion is undersized, or at least, comparable with the 
quoted statistical acceptances. Therefore, we can conclude that both of Dalitz 
decay and 7-conversion have an infinitesimal effect on ID FM, especially in the 
data retrieved by MIRACLE Lab. 



5 Results in Multi-Dimensions 



H = 1.0 



3 <ti<4 



3 <T) < 5 



2 <ti < 5 



» 6"""--"5--g----gi 

" I 



A q = 4 
□ q = 3 
• q = 2 



2 < T) < 6 



lV-"-"fr--J 



M .M 

1 <P 



10 



H = 1.0 

1 ' — i *- 

3<n<4 * 

. 8 6 8 8 88BSS" 



S t u a u g85»S.- 



4. * 



AO at BIS" 




Figure 6: FM is performed here 
in two-dimensions simultaneously and 
then compared in the different An- 
intervals. The partition process is per- 
formed for small number M v ■ M§. 
Here, FM apparently increase faster 
than the ID FM (Fig. ^ and Fig. ||) 
with the increasing of M v ■ M^, . 



Figure 7: The same as in Fig. \Q, but 
for large M v ■ M^. We notice that the 
increasing of \ogF q vs. \og(M v ■ M^) 
is upwards binding. Saturation regions 
as in Fig. |^ are observed here, espe- 
cially in the largest Arj-interval. 



For clear understanding of the dynamics of the interaction and to investigate the 
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different reasons responsible for the multiplicity fluctuations and consequently 
to explain the power-law behavior at the final state of particle production, it is 
necessary to investigate FM in multi-dimensions (0, p t , rj) |H| . Because of the 
absence of any external magnetic field during the emulsion exposure (Section 

t, it is not possible to measure pt- It remains only the two dimensions, r\ and (j> 
(1. For the multi-dimensional study of FM, we should first define the so-called 
multi- dimensional partition method. In next subsection, we introduce methods 
generally used for the two-dimensional partition. 



5.1 Methods for two-dimensional partition 

If the Ary-interval is divided into the same number of bins as A(f>, this method is 
called self-similar partition, i.e. the partition in certain dimension controls the 
method of the partition in the other one. Clearly, this method leads to total 
partition number of M 2 . A practical method to contemporaneously determine 
the partitions in two-dimensions is given by the so-called Hurst- exponents p3], 

w lnMl no\ 

Mi\ i — {1, 2} is the number of bins in i-th direction. If 7i = 1, the two variables, 
•q and <j>, will be divided in a self-similar way. It was pointed out in [^7|, |33| and 



discussed in Section 4.1 that the superposition in ^-dimension makes the effective 
partition number, ^> M v . In the contrast, the superposition in (^-dimension 
makes no change at all, i.e. M^ s — M$. These lead to Tt^ = lnM^/lnM^, 
and 



M f = e n ^ lnM *. (13) 

To study the possible abnormal behavior of FM with M, we should perform the 
analysis of FM for only one value of Ti. Otherwise, the observed trend of the 
calculated FM will always be bending upwards |28j , even if there is no fluctuation 
in the particle production. The upwards increasing of FM which is depending 
only on the partition method, can be weakened or totally removed (review Eq. 
(fl3|)), if the exponents, 7i, are given a suitable value. At this position, it is 
worthwhile to refer to the conclusion of [^9| that the exponents, 7i, applied in 
the particle production strongly depend on the energy and the mass of projectile 
and target. 



5.2 Results in two-dimensions 

F q are drawn in Fig. [(] as functions of M^M^ for the different A77 intervals. The 
partition processes have been performed in a self-similar way. The number of 
partition, M^M V , goes up to 100 for the different Arj intervals, which enable us 
to study FM for different bin sizes. Here there is no restriction on the azimuthal 
angles (A</> = 2tt and Sep e {0, 2ir}). In addition to the features concerning with 
the 2D-analysis, almost similar features as in Fig. I and Fig. § are observed 
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here, for example, the reduction of FM with the increasing of the considered Arj- 
window and the large dispersions of FM with the decreasing of A77. Obviously, 
the exponential increase of FM (Fig. [7]) can be proportionally connected with 
A77- widths (see Fig. || and Fig. ||) . The last could be understood, according to 
the facts that the increase of Ar/-width leads to an increasing of the size of the 
resulted bin Sij. These observations support the assumption that the upwards 
increasing is to be explained due to the superposition of the non -coincident sub- 



intervals of the sub pseudo-rapidity values |27|, |28| (Section [f.l| ). The partition 
numbers drawn here are smaller than the effective one, M . A ccording to Eq. 
( [l3| ) , M cS are functions of the Hurst-exponents H (Section 5A ) . Continue the 
partition process in Fig. 0, FM keep on the upwards increasing. As noticed in 
Fig. |[ saturation regions get formed in wide A?7-intervals. In bottom picture 
of Fig. the saturation curves are completely formed. 

5.3 Influences of 7 conversion on 2D FM 



As discussed in Section 4.2, most of 7-conversions are expected to take place 
within the lead foil. The measuring system MIRACLE Lab is able to perc olat e 
the measurements from most of the tracks of electron-pairs (Section || and |4.2j) . 
Generally, the contributions of 7's to multi-dimensional FM |3(], Q have to be 
taken into account. The features of MIRACLE Lab discussed above, for which 
we concluded that the effects of 7's on ID FM are to be neglected, are still 
valid also in 2D FM. Since we noticed that the interacting system, Pb+Pb, 
is intermittent in ID FM (Fig. |l| to Fig. [| and Fig. ^ to Fig. ||), we do not 
constrain to go through the estimation of 7's in 2D p5fl again. Another reason 
for this renounce is the fact that the estimation of 7's is traditionally performed 
in the following way: Adding these processes to the used simulating code and 
comparing the experimental results with the simulated ones. If the last are able 
to describe the first, the fluctuations of 7's are the responsible source for power- 
law behavior of FM j3(], [54|. In our case, It is not needed to process this method, 
since the used measuring system is obviously able to filter the measurements 
from the electron-pairs. 



6 Levy Indices and Anomalous Dimensions 

6.1 Ratios of intermittence exponents 4> q /4>2 
The relation between q- and second-order FM is given as: 

\ogF q (M {vm ) = ^ logF 2 (M {im ) + c q , (14) 

4>2 and0 g are the slopes of the relations between logF 2 (M{ I) | ( ^}) and logF q (M^ v ^) 
and the corresponding M^^y. c q are constant. From these two figures (also re- 
view Fig. [l| and Fig. || again), we notice that the exponents, 4> q , do not depend on 
^{r]\4>}- Therefore, in practice the slope ratios, 4> q /4>2, can be deduced from the 
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H- 1.0 




Figure 8: The dependence of q- on second-order FM in n-dimension is depicted. 
The lines represent the linear fitting of the experimental data. In Table \Q the fitting 
parameters are listed. We notice a positive dependence of F q on Fi for all An -intervals 
in one- as well as in two-dimensions. 
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Figure 9: The ratios of intermittence exponents, [3 q = (f> q /(f)2, given as functions of 
orders, q. The solid lines represent the fitting according to the relation fcH\). The fitting 
according to relation fiF\ ) is represented by the dotted lines. The conventional linear 
power fitting are given by the thin dash-lines. The experimental data are represented 
by the solid points. 
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Table 1: The fit parameters, 4> q /(j>2 and c q , for the straight lines drawn in Fig. [3| left 

picture for ID analysis, are listed here . 

r\ q value rj q value 

c 3 0.003 ±0.001 c 3 0.005 ±0.001 

3 /0 2 1.598 ±0.009 1-543 ±0.007 

c 4 0.006 ±0.004 c 4 0.015 ±0.003 

77 < 4 /0 2 2.354 ±0.033 < 77 < 4 /0 2 2.131 ±0.035 

c 5 0.007 ±0.014 c 5 0.031 ±0.011 

5 /0 2 3.300 ±0.102 5 /0 2 2.777 ±0.110 

c 3 0.006 ±0.001 c 3 0.001 T 0.001 

3 /0 2 1.537 ±0.004 3 /0 2 1.709 ±0.023 

c 4 0.017 ±0.002 c 4 0.004 ±0.005 

77 < 4 /02 2.083 ±0.015 2< ^ <b 4 4 /0 2 2.495 ±0.056 

c 5 0.035 ±0.004 c 5 0.011 ±0.013 

5 /0 2 2.631 ±0.039 5 /02 3.321 ±0.129 



relation log-F^Mi^i.) vs. logF 2 (M^|0i.). Fig. || depicts such relations for the 
orders, q — {3,4,5}. The lines represent the power fitting of the experimental 
data. Orders q increase, the slopes and their corresponding statistical errors get 
large. Comparisons between the fitting parameters in the different A?7-intervals 
are summarized in Table |l|. The results from the 2D-analysis are graphically 
illustrated in the right part of Fig. ^ and their corresponding parameters are 
listed in Table §. From these two tables and by taking into account the different 
statistical errors in ID and 2D, we can conclude that the corresponding slope 
ratios 0«/0 2 smoothly get larger for higher analysis dimension. In left picture 
in Fig. @, especially in the largest A?7-interval, 2 < rj < 6, we notice the exis- 
tence of two separated regions. Each of them has been fitted and depicted as 
thin line. The suddenly increasing of the ratios, q /0 2 , at the boarder between 
these two regions measures the maximum amplitude of the saturation plateau of 
second-order FM (~ 1.5). At this point, g /02 — * oo. The comparison between 
one- and two-dimensional FM requires to give 7Y-exponents, used to define the 
method of partition in two-dimensions, a suitable constant value. The value, 
H = 1, obviously cannot overcome the effects of the superposition of the sub- 
intervals (elementary collisions). Therefore, we expect that FM exponentially 
increases for Srj — > (Fig. 0). This clearly explains the non-linear behavior of 
parts of the relations shown in the right picture in Fig. ||. 

6.2 Physics of power-law behavior 

As mentioned in Section [j] and discussed elsewhere 
responsible for the power-law behavior at the final state of particle production 
can be classified into two main groups: 

1. Self-similar processes. By using Gaussian distributions p6[ |, the ratios of 
q- to second-order intermittence exponents can be expressed as: 



the sources 
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Table 2: The parameters, (j> q /<j>2 and c q , from the linear power fitting drawn in Fig. 

[jj, right picture for the 2D analysis, are given here . 

r) q value rj q value 

c 3 0.024 ±0.003 c 3 0.021 ±0.003 

3 /0 2 1.579 ±0.017 3 /0 2 1.629 ±0.023 

c 4 0.084 ±0.014 c 4 0.064 ±0.012 

^ 4 /02 2.213 ± 0.077 77 4 /0 2 2.377 ±0.083 

5 c 5 5 c 5 

05/02 05/^2 

c 3 0.021 ±0.005 c 3 0.002 t 0.006 

3 /0 2 1.592 ±0.028 3 /0 2 1.965 ±0.040 

c 4 0.072 ±0.010 c 4 0.024 ±0.012 

77 4 /02 2.226 ±0.060 ^ <7 ? <b 4 4 /0 2 3.102 ±0.081 

c 5 c 5 0.075 ±0.025 

5 /02 5 /02 4.337 ±0.170 



d 



<l 



g-S*- 1 '-!, (15) 

As discussed in |37], [38|, the implementation of Levy index, [i [^9j, leads 
to: 

d q q^-q 1 , . 

^ = 2^2-^T' (16) 

where /i has a continuous spectrum, < \i < 2, (region of Levy stability). 
Experimentally, EHS/NA22 collaboration jl(| found that the index fi, in 
principle exceeds the upper boundary, 2. As we will see later, Levy index 
within its restricted stability region evidently is not able to describe the 
experimental relations between d q jdi and the orders q. Only for 2D FM 
in the rapidity region ranging from 2 to 6, we get a stable value, /i = +0.1 
(Fig-©. _ 

Relation ( Jig ) is evidently an approximation of anomalous ratios, d q jdi. 
Obviously, this relation is not satisfactorily applicable in the distribution 
tails. Meanwhile relation (|l^) can effectively determine these tails. As 
given above, Levy index has two boundaries (drawn in Fig. [l(] as dotted 
lines): 

(a) /z = 2 leads to d q jdi = q/2 corresponding to self-similar branching 
processes represented by the well-known Gaussian distributions, and 

(b) fj, — leads d q /d2 — 1 corresponding to mono-fractal behavior re- 
sulted from second-order phase-transition from QGP to hadron-gas 
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Figure 10: Dependencies of the anomalous dimension ratios, cj> q /(j>2, with the orders, q, 
for the different An-mtervals and for the possible two dimension, r\ and <f>. The dotted 
lines represent the two boundaries of Levy stable region. The experimental results are 
represented by the open circles. The Solid lines are the solutions of Eq. |7^ ) for the 
given ji-values. 



Between these two limits, Eq. (|l5]) is not valid and should be replaced by 
the relation (|l6|). 

2. If QGP indeed is to be produced in the relativistic heavy-ion collisions, 
the interacting system has to suffer second-order phase-transition during 
its space-time evolution. In this case the observed intermittence behavior 
(in form of power-scaling behavior) can be compared with the behavior in 
a system of mono- fractal (see Fig. ^ and Fig. [H]) , 

Then we can conclude that the dependence of the ratios, 9 /</>2, on the orders, 
q, reflects essential information about the reaction dynamics In following 
sections, we will study the different intermittence ratios in relation with the 
orders, q. 
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2 <T] < 6 




Figure 11: Comparison between the second-order F M a nd the partition number, M v 



for the experimental and simulated events (Section 2.1). The disagreement between 
the data and the three models becomes larger for larger M v . In bottom part, the third- 
order FM are given as functions of the second-order ones for the same data set as in 
the top picture. The line represents the linear power fitting of the experimental data. 
Generally, we conclude that non of the three models is able to describe the experimental 
results. 



6.3 Relations between 4> q /4>2 and the orders q 

In form of a power-law, relation (^) can be re- written as: 

F q <xF%«; f3 q = cf> q /cf> 2 . (18) 

It is clear to recognize that the critical powers, f3 q , are independent on the 
partition process or on the analysis dimension. This power-law is valid in one- 
as well as in multi-dimensions [p8[ . In flt2| , the predicted second-order phase- 
transition in heavy-ion collisions has been simulated by using a specific form of 
Ginzburg-Landau model |43| |. It has been found that critical exponents, (3 q , are 
independent on the temperature, T, just below the critical one, T c . Therefore, 
the behavior of (3 q effectively characterizes the phase-transition to coherent light 
and also its aftereffects 



[3 q = (q-l) v ; v = 1.304. (19) 

In our case, the aftereffects is the ratios, 4> q /4>2, or fundamentally the fluctua- 
tions in the particle production, v are not looking like the conventional critical 
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indices describing only the behavior of quantities near the critical point. In ad- 
dition to this, they are independent on the details of the interacting system and 
therefore can effectively characterize the behavior of all measurable quantities 
|| (<V02, etc). 

Figure [| describes the powers, as functions of the orders, q. As given before, 
j3 q can directly be determined from the slopes of straight-line portions in Fig. ||. 
The slopes are given in Tables [I] and ^ for ID and 2D, respectively. The lines with 
empty squares are the fitting according to the relation Almost the same 

process, but for v = 1.0, is represented by the dash-lines. The value v = 1.0 in 
Eq. ([It]) is originally obtained from the intermittent behavior of 2D Ising-model 
[E3 and also from the thermal second-order phase-transition (Section |6.l| ) . From 
Fig. [| we also notice that the first fitting is unable to describe the experimental 
data. Although the second one results lines close to the experimental points, 
it remains unable to completely predict the experimental results. Only, 2D- 
analysis in the region, 2 < 77 < 6, produces one line coincident with the points 
representing the experimental data (solid points). 



6.4 Relations between d q /d,2 and the orders q 

Figure ^ depicts the ratios, d q /d2, in dependence with the orders, q, for the 
different A^- intervals. Th e re lations between the anomalous and Renyi dimen- 
sions are given in Section Virtually, this figure can be used to distinguish 
between the two possible sources of the intermittent behavior discussed in Sec- 



tion 3.2. As we can see, only Levy indices ranging from —1.0 to —3.5 can fit our 
experimental data. The most acceptable value, /z = +0.1, is obtained for 2D 
FM in the largest Ary-interval. The reasons and interpretations for this mono- 
fractal behavior (fi — > 0, Eq. (|l6|)) we will discuss later. Again, the predictions 
of Eq. (|l7|) are exactly the lower Levy boundary, where the line fi = is. 



7 Remarks and Final Conclusions 

The following remarks have to be considered for the evidence of the thermal 
phase-transition noticed in this article: 

1. FM have to be re-performed for higher orders and dimensions. Such a 
way, we get more lines for Fig. [9] and Fig. [lOL which might support the 
observation of thermal phase-transition. For the same reason, more mea- 
surements are also required. 

2. 2D FM have to be re-investigated by using other values of 7i. Such a way, 
we can investigate the scaling-law behavior of FM away from the non- 
canonical effects stemming only from the partition method itself (Section 

©■ 

3. A?7-windows have to re-determined to specify the regions which fulfill the 
requirement of Levy stability. The non-considered interval might include 
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leading baryons with pseudo-rapidity values close to that of the projectile 
nucleons. In that case, most of particles are of course fragmentations. 
Therefore, they behave differently than the produced ones (7r's). 

4. In forthcoming steps, more weight should be given to the smaller intervals 
and to the dynamics of the interaction. The last could be studied by us- 
ing the event-by- event fluctuations. Generally, the analysis of single events 
with huge multiplicity enable us to study different physics: geometry of in- 
teraction, space-time evolution, two-particle correlation functions, micro- 
scopic and thermodynamic collisions, QCD fluctuations, phase-transition, 
etc. @. 

Even, if we could include leading-particles, it is still early to dwell on ultimate 
conclusion. Since these non-produced particles may uncommonly comport as 
the produced ones, 7r's. Although this conclusion could, in the one hand, be 
a uproarious one, especially in light of the remarks stated before, in the other 
hand, we cannot totally neglect the demonstrative evidence noticed in Fig. [to] 
and Fig. Experimentally — in only one of Arj regions — we have a stable 
Levy index, \jl — +0.1 (Fig. || and Fig. [h]). This positive small value is standing 
as a signature for the thermal phase-transition j37j. The other values, n < 0, 
indicate that Levy indices cannot effectively distinguish between thermal and 
non-thermal branching processes. 

The stable value, /i = +0.1, is obtained from the two-dimensional analysis 
of the largest Arj region. Within this region, the effects of the superposition 
in ^-dimension are much smaller than that in other ^-intervals p9| (Section 



4.1). Within it, there are also non-produced particles, which in contrast to the 
produced particles, are not sensitive to the phase-transition. As mentioned 
above, more measurements, event- by- event analysis and higher dimensions are 
highly required to demonstrate such power-law scaling behavior and then to 
confirm the evidence of the thermal phase-transition. 
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